Josephson vortices and solitons inside pancake vortex lattice in layered 

superconductors 
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In very anisotropic layered superconductors a tilted magnetic field generates crossing vortex lat- 
tices of pancake and Josephson vortices (JVs). We study the properties of an isolated JV in the 
lattice of pancake vortices. JV induces deformations in the pancake vortex crystal, which, in turn, 
substantially modify the JV structure. The phase field of the JV is composed of two types of phase 
deformations: the regular phase and vortex phase. The phase deformations with smaller stiffness 
dominate. The contribution from the vortex phase smoothly takes over with increasing magnetic 
field. We find that the structure of the cores experiences a smooth yet qualitative evolution with 
decrease of the anisotropy. At large anisotropies pancakes have only small deformations with respect 
to position of the ideal crystal while at smaller anisotropies the pancake stacks in the central row 
smoothly transfer between the neighboring lattice positions forming a solitonlike structure. We also 
find that even at high anisotropies pancake vortices strongly pin JVs and strongly increase their 
viscous friction. 

PACS numbers: 74.25.Qt, 74.25.Op, 74.20.De, 



I. INTRODUCTION 

The vortex state in layered superconductors has a 
very rich phase diagram in the multidimensional space 
of temperature-field-anisotropy-field orientation. Espe- 
cially interesting subject is the vortex phases in lay- 
ered superconductors with very high anisotropy such 
as Bi 2 Sr2CaCu20z (BSCCO). Relatively simple vortex 
structures are formed when magnetic field is applied 
along one of the principal axes of the layered structure. 
A magnetic field applied perpendicular to the layers pen- 
etrates inside the superconductor in the form of pan- 
cake vortices (PVs) i PVs in different layers are coupled 
weakly via the Josephson and magnetic interactions and 
form aligned stacks at low fields and temperatures (PV 
stacks). These stacks arc disintegrated at the melting 
point. In another simple case of the magnetic field ap- 
plied parallel to the layers the vortex structure is com- 
pletely different. Such a field penetrates inside the super- 
conductor in the form of Josephson vortices (JVs)i2i2i4 
The JVs do not have normal cores, but have rather wide 
nonlinear cores, of the order of the Josephson length, lo- 
cated between two central layers. At small in-plane fields 
JVs form the triangular lattice, strongly stretched along 
the direction of the layers, so that JVs form stacks aligned 
along the c direction and separated by a large distance 
in the in-plane direction. 

A rich variety of vortex structures were theoretically 
predicted for the case of tilted magnetic field, such as the 
kinked lattice^^, tilted vortex chains^, coexisting lat- 
tices with different orientation^ A very special situation 
exists in highly anisotropic superconductors, in which the 
magnetic coupling between the pancake vortices in dif- 
ferent layers is stronger than the Josephson coupling. In 
such superconductors a tilted magnetic field creates a 
unique vortex state consisting of two qualitatively differ- 
ent interpenetrating sublatticesi^& This set of crossing 



lattices (or combined lattice*) contains a sublattice of 
Josephson vortices generated by the component of the 
field parallel to the layers, coexisting with a sublattice of 
stacks of pancake vortices generated by the component 
of the field perpendicular to the layers. A basic reason 
for such an exotic ground state, as opposed to a simple 
tilted vortex lattice, is that magnetic coupling energy is 
minimal when pancake stacks are perfectly aligned along 
the c axis. A homogeneous tilt of pancake lattice costs 
too much magnetic coupling energy, while formation of 
Josephson vortices only weakly disturbs the alignment of 
pancake stacks. 

Even at high anisotropies JVs and pancake stacks have 
significant attractive coupling & The strong mutual in- 
teraction between the two sublattices leads to a very 
rich phase diagram with many nontrivial lattice struc- 
tures separated by phase transitions. At sufficiently 
small c-axis fields (10-50 gauss) a phase separated state 
is formed: density of the pancake stacks located at JVs 
becomes larger than the stack density outside JVs&i^ 
This leads to formation of dense stack chains separated 
by regions of dilute triangular lattice in between (mixed 
chain+lattice state) . Such structures have been observed 
in early decoration experimentsi£*i£ and, more recently, 
by scanning Hall probed, Lorentz microscopy^, and 
magnetooptical technique>i2iS£ At very small c-axis fields 
(~ several gauss) the regions of triangular lattice vanish 
leaving only chains of stacks. 17 Moreover, there are ex- 
perimental indications^ and theoretical reasoning^ in 
favor of the phase transition from the crossing config- 
uration of pancake-stack chains and JVs into chains of 
tilted vortices. JVs also modify the interaction between 
pancake stacks leading to an attractive interaction be- 
tween the stacks at large distances^ As a consequence, 
one can expect clustering of the pancake stacks at small 
concentrations. 

Many unexpected observable effects can be naturally 
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interpreted within the crossing lattices picture. The un- 
derlying JV lattice modifies the free energy of the vortex 
crystal state. An observable consequence of this change 
is a shift of the melting temperature. Strong support for 
the crossing-lattices ground state is the linear dependence 
of the c-axis melting field on the in-plane field observed 
within a finite range of in-plane fields Ai2*i2i22i24 In an 
extended field range several melting regimes have been 
observed22i24 indicating several distinct ground states of 
vortex matter in tilted fields. Transitions between differ- 
ent ground state configurations have also been detected 
by the features in the irreversible magnetization. 25,26 

In this paper we consider in detail the properties of an 
isolated JV in the pancake lattice. We mainly focus on 
the regime of a dense pancake lattice, when many rows 
of the pancakes fit into the JV core. The pancake lattice 
forms an effective medium for JVs which determines their 
properties. The dense pancake lattice substantially mod- 
ifies the JV structure. In general, the phase field of the 
JV is built up from the continuous regular phase and the 
phase perturbations created by pancake displacements. 
Such a JV has a smaller core size and smaller energy as 
compared to the ordinary JV£ The pancake lattice also 
strongly modifies the field and current distribution far 
away from the core regional 

The key parameter which determines the structure of 
the JV core in the dense pancake lattice is the ratio 
a = A/7S, where A is the in-plane London penetration 
depth, 7 is the anisotropy ratio, and s is the period of lay- 
ered structure. The core structure experiences a smooth 
yet qualitative evolution with decrease of this parame- 
ter. When a is small (large anisotropies) pancakes have 
only small displacements with respect to positions of the 
ideal crystal and the JV core occupies several pancakes 
rows. In this situation the renormalization of the JV 
core by the pancake vortices can be described in terms of 
the continuous vortex phase which is characterized by its 
own phase stiffness (effective phase stiffness approach). 6 
At large a (small anisotropies) the core shrinks to scales 
smaller than the distance between pancake vortices. In 
this case pancake stacks in the central row form soliton- 
like structure smoothly transferring between the neigh- 
boring lattice position. 

We consider dynamic properties of JVs in the case of 
small a: the critical pinning force which sticks the JV 
to the PV lattice and the viscosity of the moving JV 
due to the traveling displacement field in the PV lattice. 
The pinning force has a nonmonotonic dependence on 
the c-axis magnetic field, B z , reaching maximum when 
roughly one pancake row fits inside the JV core region. 
At higher fields the pinning force decays exponentially 
oc exp(— yjB z /Bo). We study JV motion through the 
PV lattice and find that the lattice strongly hinders the 
mobility of JVs. 

The paper is organized as follows. Section[n]is devoted 
to the static structure of an isolated JV inside the PV 
lattice. In this section we 

• consider small c-axis fields and calculate the cross- 



ing energy of JV and PV stack Ijll B|) ; 

• consider large c-axis fields and introduce the "ef- 
fective phase stiffness" approximation, which al- 
lows for simple description of JV structure inside 
the dense PV lattice in the case of large anisotropy 
(ITTCT) : 

• investigate a large-scale behavior and JV magnetic 
field flTPjl ; 

• analyze the JV core quantitatively using numeric 
minimization of the total energy and find crossover 
from the JV core structure to the soliton core struc- 
ture with decrease of anisotropy Ijll E|) ; 

• formulate a simple model which describes the soli- 
ton core structure for small anisotropies Ijll Ffl . 



In Section IIIII we consider pinning of JV by the pancake 
lattice and calculate the field dependence of the critical 
current at which the JV detaches from the PV lattice. 
In Section ll VI we consider possible JV dynamic regimes: 
dragging the pancake lattice by JVs and motion of JVs 
through the pancake lattice. For the second case we cal- 
culate the effective JV viscosity. 



II. STRUCTURE AND ENERGY OF 
JOSEPHSON VORTEX IN PANCAKE LATTICE 

A. General relations 

Our calculations are based on the Lawrence-Doniach 
free energy in the London approximation, which depends 
on the in-plane phases </> n (r) and vector-potential A(r) 
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are the phase stiffness and the Josephson coupling en- 
ergy, A = X a b and A c are the components of the London 
penetration depth and s is the interlayer periodicity. We 
use the London gauge, divA = 0. We assume that the 
average magnetic induction B inside the superconductor 
is fixediASi The c component of the field fixes the concen- 
tration of the pancake vortices n v = _B z /$ inside one 
layer. The in-plane phases <f> n have singularities at the 
positions of pancake vortices Ri„ inside the layers, 
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The major obstacle, preventing a full analytical consid- 
eration of the problem, is the nonlinearity coming from 
the Josephson term. A useful approach for superconduc- 
tors with weak Josephson coupling is to split the phase 
and vector-potential into the vortex and regular contri- 
butions, <j) n = <j) vn + <j) rn and A = A„ + A,-. The vortex 
contributions minimize the energy for fixed positions of 
pancake vortices at Ej = and give magnetic interac- 
tion energy for the pancake vortices. One can express this 
part of energy via the vortex coordinates R« j. hi gen- 



eral, the regular contributions may include phases and 
vector-potentials of the Josephson vortices. The total 
energy naturally splits into the regular part F r , the en- 
ergy of magnetic interactions between pancakes Fm, and 
the Josephson energy Fj, which couples the regular and 
vortex degrees of freedom, 

F = F r + F M + Fj (3) 

with 
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is the magnetic interaction between pancakes^ where 

u (r, z) = exp(—z)Ei (r — z) + exp(z)Ei (r + z) , (8) 

E\(u) = (exp(— v)/v) dv is the integral exponent 
(Ei(u) « —7b — liiit + u at ji <C 1 with ~ 0.577), 
r = y/R 2 + (ns) 2 , and L is a cutoff length. 

In this paper we focus on the crystal state. If the 
pancake coordinates have only small deviations from the 
positions r|°' of the ideal triangular lattice, R„; = 
Rj ^ + u„i then F v reduces to the energy of an ideal crys- 
tal F cr plus the magnetic elastic energy Fm-eI consisting 
of the shear and compression parts, 
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(e; (e t ) is the unit vector parallel (orthogonal) to kj_), 

U t {k) = C 66 k 2 ± + C/ 44 (k), 
C/ i (k) = [/ 11 (k) + C/ 4 4(k), 



Un is the compression stiffness, Cqq is the shear modulus. 
In particular, at high fields, B z ^> $o/A 2 , we have 
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The magnetic tilt stiffness^ {7 44 (k), is given by interpo- 
lation formula, which takes into account softening due to 
pancake fluctuations, 



(k) 



C iA {k z )k\ 
2(4tt) 2 A 4 



In 1 



, (9) with r 2 



+i 



2 >, 



(10) 



A at a > A and 



r C ut ~ a/4.5 at a < A, 27 with a = \J 2/(\/3n v ) being the 
lattice constant. At finite Josephson energy minimiza- 
tion of the energy with respect to the phases at fixed 
pancake positions leads to the Josephson term in the tilt 
stiffness energy^ 

The major focus of this paper is the structure of JV 
core. This requires analysis of the pancake displacements 
and regular phase at distances r± <C A c and z <§; a, A 
from the vortex center. At these distances the main con- 
tribution to the energy is coming from the kinetic en- 
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ergy of supercurrents and magnetic screening can be ne- fied: one can neglect the field contributions in the regular 
glectcd. The structure of energy is significantly simpli- and Josephson energy terms, i.e., drop A r : 

I 
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n J 
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and use asymptotics r cut k z 3> 1 in the tilt stiffness (|10fl . 
Behavior at large distances r_i_ ~ A c and z ~ a, A is im- 
portant for accurate evaluation of the cutoff in the loga- 
rithmically diverging energy of the Josephson vortex. In 
this range the Josephson term can be linearized and one 
can use the anisotropic London theoryii 



B. Small c axis field: Crossing energy 

At small fields and high anisotropy factor pancake vor- 
tices do not influence much structure of JVs. However, 
there is a finite interaction energy between pancake stack 
and JV (crossing energy) which causes spectacular ob- 
servable effects, including formation of the mixed chain- 
lattice stateAiliiSiiSiiS We consider a JV located between 
the layers and 1 and directed along x axis with center 
at y = and a pancake stack located at x — and at 
distance y from the JV center. We will calculate struc- 
ture of the pancake stack and the crossing energy. The 
JV core structure is defined by the phases <j) n (y) obeying 
the following equation 
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with y = y/XjQ. An accurate numerical solution of this 
equation has been obtained in Ref. 0. It is described by 
the approximate interpolation formula 28 
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Interaction between the pancake stack and JV appears 
due to the pancake displacements u n under the action 
of the JV in-plane currents j n (y) (see Fig. In the 
regime of very weak interlayer coupling the energy of the 
deformed pancake stack is given by 



V— j n (y)u n (15) 
^— ' c 




FIG. 1: Configuration of the pancake stack crossing the 
Josephson vortex: (a) the stack located in the center of JV 
core and (b) the stack located at a finite distance y from JV 
center. 



where U M {k z ) = 2(4 ^V ln (} + T^qbg 
netic tilt stiffness of the pancake stack, 



is the mag- 
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with p n (y) = d(j) n (y) / dy being the reduced superfluid mo- 
mentum, and the JV phase cf> n (y) is given by approximate 
formula l|14fl . In particular, p n (0) = —C n /(n— 1/2) with 
C„ — > 1 at large n. Using the precise numerical phases 
<t>n{y) w e obtain the interpolation formula 

C n « 1 - 0.265/((n - 0.835) 2 + 0.566), 

giving Ci 0.55 and C2 ~ 0.86. At large distances from 
the core, n, y ^> 1, p n {y) is given by 
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Displacements in the core region have typical wave 
vectors k z ~ ir/s. In this range one can neglect /in- 
dependence of UM{k z ), 
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and rewrite Eq. (|15|l as 
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We neglected weak logarithmic dependence of the tilt 
stiffness on displacements and the parameter r w in Um is 
just a typical value of \u n +i — u n \. Minimizing this energy 
with respect to u n , we obtain the pancake displacements 
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with v n (y) ~ \u n (y) - u n -i(y)\ and crossing energy at fi- 
nite distance y between the crossing point and the center 
of JV core 
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with 



A x (V) = E [ Pn ^)] 2 = E ^ ~ C0S (<?Wl ~ 0n)) : 
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where the second identity can be derived from Eq. (|13|l . 
In particular, for the pancake stack located at the JV 
center, A x (0) = 2 (exact value) and 
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The maximum pancake displacement in the core is given 
by 
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At large distances, y 7s, using asymptotics A x (y) 
Tt/Ay, we obtain 
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Using numerical calculations, we also obtain the follow- 
ing approximate interpolation formula for the function 
A x (y), valid for the whole range of y, 
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Eqs. (H3> and JS) dete rmine the crossing energy at finite 
distance y between the crossing point and the center of 



the JV core. Below we will use this result to calculate the 
pinning force which binds the JV to the dilute pancake 
lattice. 

C. Large c axis field. Approximation of the 
effective phase stiffness 

At high c axis fields pancakes substantially modify 
the JV structure. Precise analysis of the JV core in 
the pancake lattice for the general case requires tedious 
consideration of many energy contributions (see Section 
III El below) . The situation simplifies considerably in the 
regime of very high anisotropy 7 3> A/s. In this case 
one can conveniently describe the JV structure in terms 
of the effective phase stiffness, which allows us to re- 
duce the problem of a JV in the pancake lattice to the 
problem of an ordinary JV at B z = 0. In Ref. this 
approach has been used to derive JV structures at high 
fields B z 3> $ /47rA 2 . The approach is based on the 
observation that smooth transverse lattice deformations 
u t „(r) produce large-scale phase variations <j) vn {r) with 
V0. url = 2im v e z x u tn . This allows us to express the 
transverse part of the elastic energy, F v -t, in terms of 
<?Wr): 
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with the effective phase stiffness J v (B z ,k), 
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Replacing the discrete lattice displacements by the 
smooth phase distribution is justified at fields B z > 
$o/(7 s ) 2 - The structure of the JV core is determined 
by phase deformations with the typical wave vectors 
k±_ ~ I/7S < 2/ A and k z ~ ir/s > l/r w . In this range 
the vortex phase stiffness is k-independent, similar to the 
usual phase stiffness, 
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The phase stiffness energy (|2*2*|l has to be supplemented 
by the Josephson energy. In the core region we can ne- 
glect the vector-potentials and write the total energy in 
terms of 6 r7 , and d>„ as 
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where, again, <p n = <p rn + 4> vn is the total phase. 

We now investigate the core structure on the basis of 
the energy jS5)l. Eliminating the regular phase, (p rn = 
4>n ~ i'vni an d varying the energy with respect to <p vn at 
fixed 4> n , we obtain the equation 



which gives 



JA(f> n - (J v + J) A(f> vn = 0, 
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Substituting this relation back into the energy 1)25)1 . we 
express it in terms of the total phase 



with the effective phase stiffness J e ff, 
J~s = J -1 + J" 1 or J cff 
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Note that the smallest stiffness from J and J u dominates 
in Jcfff^i From Eq. (|27|l we obtain equation for the equi- 
librium phase 

JeffV^n + [sin (0 n+1 - i/» n ) - sin - n _i)] = 0, 

(29) 

which has the same form as at zero c axis field, except 
that the bare phase stiffness is replaced with the effective 
phase stiffness J e g. For the Josephson vortex located be- 
tween the layers and 1 the phase satisfies the conditions 
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(30) 



Far away from the nonlinear core the phase has the usual 
form for the vortex in anisotropic superconductor 



4>n{y) ~ arctan 
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where the effective Josephson length 

A,/ = y/J s/Ej = Ajo/v/1 + B z /B\ 

determines the size of the nonlinear core. Therefore, at 
low temperatures the JV core shrinks in the presence of 
the c axis magnetic field due to softening of the in-plane 
phase deformations. A number of pancake rows within 
the JV core can be estimated as 



A jo / ln(r cut /r w ) 
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At B z > B\ it is almost independent on the field. An 
approximate solution of Eq. (|29|l is given by Eq. 1)14)1 . 
where the bare Josephson length Ajo has to be replaced 
by the renormalized length A,/. 



The JV energy per unit length, £jv, is given by 

L 



£.jv = iry/EjJ cfi In 



(33) 



where L is the cutoff length, which is determined by 
screening at large distances and will be considered be- 
low, in Section Iff Dl From Eqs. iffify and l|2*H)l we obtain 
that the partial contribution of the vortex phase in the 
total phase, 



B, 



B z + B x 
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continuously grows from at B z <C B\ to 1 at B z ^S> 
B\. From the last equation one can estimate pancake 
displacements 
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The maximum displacement in the core can be estimated 
as 

2.2A 2 
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At B z 3> B\ this equation can be rewritten in the form 
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which shows that condition for applicability of the linear 
elasticity, u X fi(0) < 0.2a, is satisfied if 7 > 3A/s. Eqs. 
(131)1 and 1)33)1 describe smooth evolution of the JV struc- 
ture with increase of concentration of pancakes starting 
from the usual vortex at B z = 0. It is quantitatively 
valid only at very high anisotropics 7 3> A/s and at 
low temperatures, when one can neglect fluctuation sup- 
pression of the Josephson energy. Thermal motion of 
the PVs at finite temperatures induces the fluctuating 
phase 0n,n+i and suppresses the effective Josephson en- 
ergy, Ej — * CEj where C = (cos (f> n _ n+ i^ . This leads to 

reduction of the JV energy and thermal expansion of its 
core. 

In the range <&o/(7s) 2 < B z < B\ the "crossing en- 
ergy" regime of Section III Bl overlaps with the applica- 
bility range of the effective phase stiffness approxima- 
tion. To check the consistency of these approximations 
we calculate correction to the JV energy at small fields 
summing up the crossing energies and compare the result 
with prediction of the "effective phase stiffness" approx- 
imation. The correction to the JV energy is given by 
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where L y = Mb is the long-range cutoff length. On the 
other hand, Eqs. J2EJ) and give at B z -C B x 



5£. 



B, 



JV 



(39) 



This result is identical to Eq. (|38[) except for expressions 
under the logarithms, which are approximate in both 
cases. 



D. Large-scale behavior. Screening lengths 

In this Section we consider the JV structure at large 
distances from the core, n > 1, y > A,;. At large 
distance screening of supercurrents becomes important 
and one can not neglect the vector-potential any more. 
At these scales the phase changes slowly from layer to 
layer so that one can expand the Josephson energy in 
Eq. with respect to the phase difference and use 
the continuous approximation, n +i — 4> 



2tts a 
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sE, 



2tt 



-A: 



(40) 

This reduces the Lawrence-Doniach model defined by 
Eqs. J3J|, JSJ, and © to the anisotropic London 



model. Within this model the JV structure outside the 
core region has been investigated in detail by Savel'ev et. 
al4i. In this section we reproduce JV structure at large 
distances using the effective phase stiffness approach. For 
the vortex energy one still can use Eq. I|22[) with the full 
fc-dependent phase stiffness (12. '{t . Within these approxi- 
mations the energy © is replaced by 

F[<j> r ,<j> v , A] = F r [(f> r) A x ] + F v - t [<l>v] + Fj[<j) r + <j) v ,A z ]. 

(41) 

Varying the energy with respect to A, we obtain 



^Vx^- A ± + A 2 V 2 A ± = 0, 
P-V z cf>-A z + \ 2 C V 2 A Z = 0. 

ZTT 



(42a) 
(42b) 



It is convenient to perform the analysis of the large-scale 
behavior in k-space. Solving linear equations (|42a|l and 
(I42b|l using Fourier transform, 



A ± = 

A z = 



2vr 1 + A 2 fc 2 ' 
2tt 1 + A 2 fc 2 ' 



(43a) 
(43b) 



and excluding A, we express the energy in terms of 
phases 



F 



rf 3 k 

(2tt) 3 



J A 2 fc 2 , J„(k) a sE., Agfc 2 2 



Following the procedure of the previous Section, we eliminate <f> r , minimize the energy with respect to <f> v , and obtain 
the energy in terms of the total phase 



F = 



d 3 k 



Jcff(k) 2 , sEjjk) 2 

(V_L<P) + r (V z 0) 



2.s 



(44) 



where the effective phase stiffness, J e gp(k), and the effective Josephson energy are given by 



EjQl) = E 



A 2 fc 2 
Agfc 2 
7 1 + A 2 fc 2 



(45a) 
(45b) 



In the case of the JV, minimization with respect to the phases has to be done with the topological constrain, 
V z V y — V y V 2 = 2TrS(y)S(z), which gives 



27rifcj / J cff (k) 



2mk y (1 + A 2 fc 2 ) 



J cff (k)fc 2 + s 2 Ej{k)k 2 k 2 (1 + X 2 c k 2 + X 2 k 2 {\ + w(k)) ' 
2irik z s 2 Ej(k) _ 2nik z (l + A 2 fc 2 ) 



J off (k)fc 2 + s 2 Ej(k)k 2 fc 2 (1 + A 2 fc 2 + A 2 fc 2 (l + u,(k))) ' 



(46a) 
(46b) 
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and 



-JV 



1 

2 

J 

2s 



d 2 k 



S Ej(k)J cS (k) 
d 2 k 



\- 2 +j 2 k 2 + k 2 {l + w(k)) 



with 



;(k) = J/J (k) = 



2/i 



X 2 k 2 /2 + ln(l + k 2 z r 2 cut ) 



(47) 



(48) 



and ft, = 47m„A 2 . The integration has to be cut at k z ~ 7r/s. In addition, integration with respect to k y is typically 
determined by k y ~ fc z /7 so that one can neglect in tu(k) the term \ 2 k 2 /2 ~ \ 2 k 2 /2^ 2 1 coming from the shear energy, 
in comparison with the tilt energy term In (l + k 2 r 2 Mt ) and the JV energy reduces to 

2hk 2 



£jV ~Ts 



U2 r 2 > 
"a ' cut/ 



-1 



7TJ 



dk. 



13 Jo ^A- 2 + fc 2 (1 + 2h/ In (1 + rlJ (fc z " 2 + r 2 ))) ' 



(49) 



r 



A similar formula has been derived in Rcf. [U This for- 
mula shows that the small-fc z logarithmic divergence in 
the integral cuts off at k z = max(A _1 , a -1 ). To repro- 
duce JV energy at B z — the upper cutoff has to be 
chosen as k c — 2.36/s. 

Lets consider in more detail the case of a large c axis 
field h 3> 1 , where JV structure is strongly renormalized 
by the dense pancake lattice. Formula for the JV energy 
simplifies in this limit to 



To estimate this integral, we split the integration region 
into two intervals, l/r w < k z < n/s and Tr/a < k z < 
l/r w , and obtain 



-JV' 



7TJ 



In 



0.2a 



3/2\ 



£ 



ttJ 



7sV2/i 



Note that the long-range contribution to the energy 
scales as a logarithm to the power 3/2. 



J 



1. Magnetic field of Josephson vortex 
Using Eqs. (|46|) and l|43|) we obtain for the JV magnetic field (see also Ref. ITU 

B x (k) = 



(50) 



I + A2fc2+A 2 fc2(l +U ;(k))' 

where w(k) is given by Eq. (|48|l . Let us consider the case of large magnetic fields B > B\. In a wide region 
A < \J y 2 + ('jz) 2 < A c , the magnetic field in real space is approximately given by 

2X 2 hk 2 ' ' 



f dk dk 

B x (y, z) « $ / ^ 2 Z exp (£fc y j/ + ifc 2 z) ( A 2 A; 2 + 



ln(l + a 2 fc 2 /20) 



(51) 



27rAA c V2ft Jo 



afc 2 — : exp 



2hk z \y\ 



1y/\a. (1 + a 2 k 2 /20) 



cos (k z z) 



One can estimate from this expression the JV maximum and this field decays at the scale ~ a/4.5 in the z direc- 
field as tion and at the scale 7a 2 /20A in the y direction. The 

3/2 



fi,(0,0)« - *° 7= flnf) " (52) 
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magnetic flux concentrated at this region is estimated 
as $ w $ /(l + 2.8/i 2 ). The residual flux, $ - is 
distributed over the pancake lattice at much larger dis- 
tances. 

Due to the elasticity of the pancake lattice, the be- 
havior at large distances is very unusual. The limiting 
expression for B x (k) at k — > is given by 



S x (k)«$ 1 



hkl 



kl/A + kll{2.Eh) 



(53) 



Formally, the total flux of JV is given by the limit 
$ = lim BJk) 

However this limit depends on the order of limits limfc^o 
and lim fc ^ 

lim lim Bj,k) = — - w — 

fe.-o^-o v ' 1 + 2.8/1 2 2.8/i 2 

lim lim BJk) = <t> 

fc„— »0fe*->0 

This apparent paradox can be resolved by calculating the 
field distribution in the real space 



1 - 



1 



$o<5(r) (54) 
$0 4y 2 - z 2 /h 



This expression clearly shows that the screening is in- 
complete: the field at large scales has a slowly decaying 
1/r 2 tail. The magnetic flux through the large size box 
L y x L z is given by 

*(L Jf> L,)=*o (l-(l-^)2arcta»(*^ 

and the limiting value of the total flux at L y , L z — > 00 
depends on the aspect ratio L y /L z . 



E. Quantitative analysis of the core structure 




i=0 j=i 

FIG. 2: Displacements of the pancake rows in the JV core. 



of the displacements of the vortex rows u n ^. Different 
representations for the magnetic interaction between the 
vortex rows UMr(u n ,i — u m j,n — m) are considered in 
Appendix [X] We will operate with the phase pertur- 
bation </> n (r) with respect to equilibrium phase distribu- 
tion of the perfectly aligned pancake crystal. We split 
this phase into the contribution, averaged over the JV 
direction (x axis), 4> n (y)i and the oscillating in the x 
direction contribution, <p n (x,y). Pancake displacements 
induce jumps of the average phase at the coordinates of 
the vortex rows Yi, 4> n (Yi + 0) — 4> n {Yi — 0) = 2iru n ^/a. 
The oscillating phase induced by the row displacements 
becomes negligible already at the neighboring row. This 
allows us to separate the local contribution to the Joseph- 
son energy coming from <j) n (x,y) (see Appendix |B|) and 
reduce initially three-dimensional problem to the two- 
dimensional problem of finding the average phase and 
row displacements. Further on we operate only with the 
averaged phase and skip the accent "— " in the notation 
4>n{y)- We again split the total phase into the continu- 
ous regular phase, 4> rn (y), and the vortex phase, 4> vn (y), 
4>n(y) = 4>m{y) + 4>vn{y\ The vortex phase is com- 

posed of jumps at the row positions Yi, 



The simple "effective phase stiffness" approximation, 
described in Section III CI is only valid if 7 is significantly 
larger than A/s. In BSCCO, even at low temperatures, 
7 is at most 3 — 4 times larger than A/s. Moreover, it 
always approaches A/s at T — > T c . In this section we 
extend our analysis to the region 7 ~ A/s. We consider 
JV structure at low temperatures and not very small c 
axis field, B z > $o/(7s) 2 . The structure of JV core is 
completely determined by the displacements of pancake 
vortices and phase distribution. The equilibrium pancake 
displacements depend only on the layer index and on co- 
ordinate, perpendicular to the direction of the vortex (see 
Fig.Ejl. Therefore, the energy can be expressed in terms 



<pvn(y, u n ,i) = y~] u„,iO (Yi - y) , (55) 

n. L — » 



where (y) is the step-function (6 (y) = 1 (0) at y > 
(y < 0)). In the effective phase stiffness approach of Sec- 
tion lll Cl we used a coarse-grained continuous approxima- 
tion for this phase. We neglect a;-dependent contribution 
in the regular phase, which is small at B > $o/(7s) 2 - 
Collecting relevant energy contributions, we now write 
the energy per unit length in terms of the regular phase, 
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i(y), and the row displacements, u„ t i, 

J / d<j> rn \ 
2 \dy~j 

- U~Mr{u n ,i - ) a - to) 



+ Ej (1 - cos 0„ + i - <j) n )) 



(56) 



• £ Josc (u n+1 ^-u n . il <p n +i,i-<f>n,i) is the localJoseph- 
son energy due to the oscillating component of the 
phase difference (see Appendix [5J) with 



(r^ + TT^-^^^-e^-y) 



being the external phase at i-th rows and n-th layer. 



where The energy l|56|l describes JV structure at distances r± <C 

j. i \ j. i \ , j. i \ • x-u i. ^ i i. A c and z <C a, A from its center. 

• 9n\y) = (f>rn{y) + PvnW, u n ,i) is the total phase; 

To facilitate calculations we introduce the reduced co- 

• U Mr {x n.i 2Cm,j-) 

Yi t j , n—m) is the magnetic interac- ordinates 
tion between the vortex rows separated by distance 
Yij = Yi — Yj = b(i — j) (see Appendix |A"|) . 



UMr(u n ,i ~~ u m,j, ■ ■ ■) — UMr(Xi-j + u n,i ~ u m,j, ■ ■ ■) 
— UMr(Xi-j, . . .). 

is the variation of this interaction caused by pan- 
cake row displacements, Xi = for even i and 
Xi = a/2 for odd i. 

UMr(x,y,n) is periodic with respect to x, 
U Mr (x + a, y, n) = U Mr (x, y, n); 

I 



y 

y = — , v r , 

7s 



and represent the energies in the scaling form. Magnetic 
interaction between the rows we represent as 



U M r{x,y,n) = -r^VMr -, -,n), 
A z a a 



where 



V Mr (x,y,n) = ( 5n - -r- exp ( -^i ) ) In [1 - 2 cos 2ttx exp (-2n\y\) +exp(-47r|y 



2A 



sA x -> 

0^2 



\/y 2 + (x — to) 2 s\n\ 



2a 2 ^ ~\ A/a ' A 



and w(r, z) is defined by Eq. (JHJ). Note that at x,y — > VMr(x,y,n) remains finite for b / because, logarithmic 
divergency in the first term is compensated by logarithmic divergency of the to = term in the sum. At p 2 = 

x 2 + y 2 ^ 0, using asymptotics it(p, z) exp(— z) ^— 7e — In (fz)) + exp(z)_Ei (2z) with je = 0.5772, we obtain the 

limiting value of Vmt{x, y, n) at n 



V A/r (0,0,7i) = ^ 



exp 



A 



In 



lTr 2 s\n\X 



7_e + exp 



s\n\ 



Er 



2s\n\ 



111 — 1 



am s\n\ 



The local Josephson energy can be represented as 

£josc{u, 4>) = 2Ejacos{cj))J(u/a), 

where J(v) is dimensionless function, J(v) « jv 2 In at v -C 1 (see Appendix |B|) . In reduced units the total 
energy takes the form 



£j/ejo = ^2 J d y 



1 fd<P 
2 



dy 



n^fsa 
~2A 2_ 



1 — cos ( 



^ J(Vn+l,i - V n ,i) C0S(<^„_|_1 i - 4> n i ) 



(57) 
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where ejo = Ejjs = J/-fs is the JV energy scale, (j> n {y) 
energy, we obtain equations for v n i and <j) n 



i(y) - 2nJ2i «n,i© (Yi - y) and a 7 = 0/73. Varying the 



+ ^TI ^2 FMr(v n ,i ~ "mj, — ,Tl ~ 77l) + V] COS(4> n ,i - (j>n+S,i)Fj(v n ,i ~ V n+S ,i) = 0, (58a) 



y^rn 



dy 



r,2 



^n+l 



5=±1 



,) - sin (</> n - n _i) = 0. 



(58b) 



Here J^Mrix, y, n) = — V x Vmt{x, y, n) is the magnetic interaction force between the vortex rows 

2tt\ 2 /„ s ( s\n\\\ sin27rx 



^Mr(x,y,n) 



s / sn 
5n -2X eXP (- — 



sA \ -t x — to 

°* m (x-m) 2 + y 2 



exp 



cosh27r|y| — cos27rx 
\J(x — mf +y 2 + z 2 



A/a 



and T,j(v, (j>) — —dj(v)/dv w — (%/2)v ln(0.235/w) at u <C 1. The derivative of the regular phase has jumps at the 
positions of the rows 



2a 



™ (Yj +0) 3?C^i - 0) = — V" i7(«n+<5,i - v ni ) sin((j) n+ s,i - 4> n ,i). 

n dy ^ ^— ' 



7 s 



(59) 



<5=±1 

r 



To find JV structure at low temperatures one has to solve ready a very weak influence on JV properties. The sum 



Eqs. |(5HaJ|, H58bf) . and with condition |gDJ. 

Let us consider in more detail magnetic interactions be- 
tween vortex rows, i.e., the term with J- Mr m Eq. (|58a|l . 
Firstly, one can observe that the dominating contribu- 
tions to the sum over the layer index to and row index j 
come from rows in the same layer, m = n, and rows in 



over the layers in / t nt [v n ,i] is determined by large num- 
ber of the layers of the order of a/s or A/s. If we consider 
layer n close to the JV core, then the interaction force 
with row in the layer m, in the same stack with n <C 
m < a/s, A/s is given by Tmt (« n ,t - «Vi,0,n- to) w 
— {vn,i — "m.i) /2(m — n). Interactions with remote lay- 



the same of stacks, j = i. The former sum determines ers give large contributions even if displacements in these 
the shear stiffness, while the latter one determines the layers are small, v m ^ <C v n ,i- A useful trick to treat this 



magnetic tilt stiffness. 



^ ^Mr (Vn,i ~ «m,j i ~^"' n_TO ) ~ /shear ["n,i] + /tilt [«n,i] 

(60) 



with 



/shear [^n,i] — ^ ^ A/r {Vn.i Vn,ji ) 0) 

/tilt fara.i] 



^ FMr{v n ,i - u m ,i,0,n- to) 



(61) 
(62) 



situation is to separate interaction of a given pancake row 
with the aligned stack pancake rows: 

/tilt [v n ,i] = Tmt (Vn,i ~ V m ,i, 0, n - to) + / cagc (l>„,j) 

with 

Fhlr (v n ,i - "m,., 0, 71 - to) 

= T M r {v n ,i - v mti ,0, n-m) - T M r {v n ,i, 0, n - to) 



and 



/cage («) = X! ( W ' ' TO ) 

m/0 



The sum over the rows in /shear [i>n,i] converges very fast 
and effectively is determined by the first two neighbor- 
ing rows. Note that skipping the terms with to =/= n in 
/shear [vn,i] is completely justified in the limit a < A but is the interaction force of a chosen pancake row with the 
leads to overestimation of the shear energy in the limit aligned stack of rows ( "cage" force) , for which we derive 
a > A. However in this limit the shear energy has al- a useful representation 

I 



/cage (V) = ~ ^ ■ 



47T sin (2ttIv) 



i=i y] (a/A) 2 + (2vr/) 2 ( y] (a/A) 2 + (2irl) 2 + 2nl 
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Atu^;l^;A/a this equation gives / ca gc (v) ~ —v ln(0.433/u). J- Mr {v n .i — v m ,i, 0, n — m) behaves as v m ^/2(m — n) 
at large m and decays with increase of m much faster than Tmt i v n,i ~ w m) i,0,n.— m). The same splitting can be 
made in the magnetic coupling energy: 



- ^ VMr{v n - V m , n-m) = ^ ( v Mr{v n -v m ,n-m) - V Mr (v n , n - m) - V Mr (v m ,n - m)) 

\m\ > \n\ 
+ ^ ^cagc (^n) 



2 

m=jin |m|>|n| 



with 



«cagc (v) = ^ VMr{v, n) 

2 (1 - cos(2ttZu)) 



E 



J=i l\l{a/\) 2 + (2nl) 2 (j(a/\) 2 + (27rlf + 2nl 



This function has simple asymptotics atv^;l^;A/a, 
«cagc(«) « (« 2 /2)ln(0.713/«). 

We demonstrate now that in the limit 7 » A/s Eqs. 
(I58a[l and (|58b|l reproduce the JV structure obtained 
within the effective phase stiffness approximation. One 
can show that in this limit the local Josephson energy 
influences weakly the JV structure. We calculate correc- 
tion to the JV energy due to this term in the AppcndixlUl 
The dominating contribution to the magnetic interaction 
between the pancake rows, Y, m ,j ^Mr{v n a-v m ,j, ^^i n ~ 
m), comes from the tilt force (|62() . which with good ac- 
curacy can be described by the cage force / ca gc ( v n,i) 
in the limit of small i> n ,i, /cage (v) ~ —v\n(C/v) with 
C « 0.433. Further estimate shows, that the term 
Hm^n Fmt (vna - v m ,i, 0, n - m) in Eq. amounts to 
the replacement of numerical constant C under the loga- 
rithm by slowly changing function of the order of unity. 
Because of slow space variations, the discrete row dis- 
placements v n< i can be replaced by the continuous dis- 
placement field v n (y). Within these approximations Eq. 
(|58a|l reduces to 

_ , 7sa C 
V^„«^ n ( y )ln-^ 

Replacing v{y) by the vortex phase 4> vn (y) obtained by 
coarse-graining of Eq. JSSJ), v n (y) = (b/2-K^s) V y 4> vn (y), 
we obtain 

Vy(Prn ~ V 'y(f)vn 

Note that we replaced v n (y) in lri[C/v n (y)} by its typ- 
ical value and absorbed the logarithmic factor into the 
definition of B\ (|24() . From the last equation we obtain 

4>m = {B\/B z )(f> m axid 




Therefore Eq. I|58bjl reduces to 

-- — 1 —r-~TT' + sin 0An ~ <An+l) + sin (<f> n - <j> n -l) = 0, 

1 + bZ dy 

which is just dimensionless version of Eq. 129(1 . 



1. Numerical calculations of J V core structure. Crossover 
to solitonlike cores. 



To explore the JV core structure, we solved Eqs. I|58af) 
and (|58b|l numerically for different ratios A/7S and dif- 
ferent magnetic fields. We used a relaxation technique 
to find the equilibrium displacements of the pancake 
rows and the continuous regular phase. Typically, we 
solved equations for 20 layers and the in-plane region 
< y < 20. 

To test the "effective phase stiffness" model and to cal- 
culate uncertain numerical factors, we start from the case 
of large anisotropies, 7 X/s. Figure [3] shows the grey 
level plots of the cosine of the phase difference between 
two central layers of JV, cos 9, 9 = <j>i — (f>o = 2<f>i, for 
A = 0.27s. Fig.^Jshows y dependence of the total phase 
difference 9 and the contribution to this phase coming 
from the regular phase for the same parameters. As one 
can see, at B > <f>o/(7s) 2 the core region covers sev- 
eral pancake rows. At high fields the core size shrinks 
so that the number of rows in the core does not change, 
in agreement with the "effective phase stiffness" model. 
From Fig. 0] one can see that the fraction of the regular 
phase in the total phase progressively decreases with in- 
crease of magnetic field. For the field 8<&o/(7s) 2 we also 
plotted &(y) dependence from the "effective phase stiff- 
ness" model, assuming B\ = 2.1$ /47rA 2 . One can see 
that the numerically calculated dependence is reasonably 
well described by this model. 
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B [in units O /(ys) 2 ] 



8 




FIG. 3: Grey level plots of cosine of the phase difference be- 
tween two central layers of JV, cosO, for A = 0.27s and sev- 
eral magnetic fields (dark regions correspond to cosO ~ — 1 
and white regions correspond to cosO ~ 1). The total size of 
displayed region in the horizontal direction is 67s. One can 
see that in this regime several pancake rows fit inside the core 
region. At high fields the size of the core shrinks so that the 
number of pancake rows inside the core remains constant. 



We now extend study of the core structure to moder- 
ate anisotropics ~ A/s. At Fig. we plot the maximum 
pancake displacement it max in the core region normalized 
to the lattice constant a as function of magnetic field 
for different A/7S. The maximum displacement approx- 
imately saturates at a finite fraction of lattice constant 
at high field (at high A/7S one can actually observe a 
slight decrease of u max /a with field). Fig.||)]shows depen- 
dence of u m ax/a on A/7S for fixed field B = 10&o/ (js) 2 . 
Dashed line shows prediction of the "effective phase stiff- 
ness" model given by Eq. (|37Jl . One can see that this 
equation correctly predicts maximum displacement for 
A/7S < 0.35. Important qualitative change occurs at 
A/7S > 0.35, where the maximum displacement it m ax(0) 
exceeds a/4. This means that the pancakes initially be- 



© 




© 



© 



© 




B = 8* /(Tfsr 



- Effective phase stiffness | 



FIG. 4: Coordinates dependence of the phase difference be- 
tween two central layers, = <j>\ — <j>o, for A/7S = 0.2 and 
different magnetic fields. Circles connected by dotted lines 
represent total phase difference, solid lines show contributions 
from the regular phase. Jumps of the total phase difference 
at the positions of pancake rows are caused by pancake dis- 
placement and represent the vortex phases. In the lower plot 
dashed line represents prediction of the "effective phase stiff- 
ness" model with B A = 2.1$ /47rA 2 . 



Maximum pancake displacement in the core 




B (in units <S>q/4tzX ) 



FIG. 5: Field dependencies of the maximum pancake dis- 
placement in the core at different ratios A/7S (the curves are 
labelled by this ratio) . 
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FIG. 6: Dependence of the maximum pancake displacement 
in the core on the ratio A/7S at B = W^o/AivX 2 . Dashed line 
represents prediction of the "effective phase stiffness" model 



longing to the neighboring stacks become closer than 
the pancakes belonging to the same stack. This can be 
viewed as switching of the vortex lines in the central layer 
of JV. This switching is clearly observed in Fig.0 which 
shows pancake displacements in the central row of pan- 
cake stacks and its neighboring row for two values of the 
ratio A/7S, 0.3 and 0.5, and several fields. For X/js = 0.5 
configuration of the pancake rows in the central stack is 
very similar to the classical soliton ("kink") of the sta- 
tionary sine-Gordon equation: the stacks smoothly trans- 
fer between the two ideal lattice position in the region of 
the core. Simplified approximate description of such so- 
litionlike structure in the case 7s < A is presented below 
in Sec. Ill Fl Fig. [S] shows distribution of cosine of inter- 
layer phase difference between two central layers, cosO. 
As one can see, at A/7S = 0.3 there are still extended 
regions of large phase mismatch in the JV core (dark 
regions), while for A/7S = 0.5 these regions are almost 
eliminated by large pancake displacements in the core. 



configuration a strong phase mismatch between the two 
central layers is eliminated, which saves the Josephson 
energy in the core region. On the other hand, large pan- 
cake displacements lead to greater loss of the magnetic 
coupling energy. To describe this soliton structure, we 
consider a simplified model, in which we (i) keep only 
displacements in the central row u„.o = v m (ii) use the 
cage approximation for magnetic interactions and (iii) 
neglect the shear energy. All these approximations are 
valid close to the JV center. At n > we redefine dis- 
placements as v n — > l + v n . The redefined displacements 
depend smoothly on the layer index, - u n C 1, so 
that one can replace the layer index n by continuous vari- 
able z = ns, u n — > u(z), and use elastic approximation 
for the Josephson tilt energy: 



dz 



nEjsCj ( du\ nJa 



2a \ dz 



Wcage (J 



(63) 



where the logarithmic factor Cj is estimated as Cj 



In 0.39(a/s) \du/dz\ and 



^ 1 — cos (2nlv 

Wcagc(w) W 2^ (27r) 2 Z 3 



1=1 



is the magnetic cage. For estimates, we replace Cj by a 
constant substituting a typical value for du/dz under the 
logarithm. The equilibrium reduced displacement v = 
u/a is determined by equation 



(64) 



Its solution is implicitly determined by the integral rela- 
tion 

dv 



A^z 

/2 y/v cagc (v) A 



F. Simple model for solitonlike cores at moderate 
anisotropies 

In this Section we consider the structure of the JV 
core for moderate anisotropies 7s < 0.5A and high fields 
B > <I>o/47rA 2 . Estimate (|37l) and numerical calculations 
show that at sufficiently small anisotropy, 7 < 0.5A/s, 
the maximum displacement in the core region exceeds a 
quarter of the lattice spacing. This means that distance 
between displaced pancakes belonging to the same vortex 
line, 2u x fi(0), becomes larger than the distance between 
pancakes initially belonging to the neighboring lines, 
a — 2u x fi(0) . This can be viewed as switching of the vor- 
tex lines in the central layer of JV. At lower anisotropics 
pancake stacks in the central row acquire structure, simi- 
lar to the soliton of the stationary sine-Gordon model: in 
the core region they displace smoothly between the two 
ideal lattice position (see Fig.[7]for A/7S = 0.5). In such 



with A = J£. Therefore a typical size of the soliton is 
given by 

Zs ~ A/7 

and the applicability condition of this approach z s 3> s 
is equivalent to 7s A. In fact, an accurate numerical 
calculations of the previous section show that the core 
acquires the solitonlike structure already at 7s < 2A. 
The core energy is given by 



any/2JEj£j 



A 



- / dv^Jv cage (v) 



At a <C A numerical evaluation of the integral gives 
Jq dvy/ v cagc (v) « 1.018/27T and we obtain 



J 2 A 
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FIG. 7: Structure of the pancake-stacks row in the center of JV (big circles) and its neighboring row (small circles) for A/7S = 0.3 
and 0.5 and several values of the magnetic field. At A/7S = 0.5 pancakes in the central row form lines smoothly transferring 
between two ideal lattice position (solitonlike structure) . 



B (in units ^/(ys) 2 ) 
My$ 1 2 4 




FIG. 8: Gray-level plots of cosine of interlayer phase difference between the two central layers of JV for A/7S = 0.3 and 0.5 and 
different fields. For A/7S = 0.3 the JV core covers roughly three pancake rows, while for A/7S = 0.5 it shrinks to one pancake 
rows. In the second case the regions of suppressed Josephson energy are practically eliminated. 



We also estimate the shear contribution to the periodic 
potential v cagc (v) 

2A 2 2 (1 - cos 2irv) exp (-y/3ir) 
a (l + cxp(- V / 37r)) 
A 2 

« 0.017 — (l-cos27ro). 
cr 

It occurs to be numerically small and only has to be taken 
into account when a becomes significantly smaller than 



A. 



It is important to note that the described model does 
not provide precise soliton structure. At distances z > 
A/7 displacements in other pancake rows become compa- 
rable with displacements in the central row. In fact, at 
distances z » A/7 the displacements should cross over 
to the regime, described by the "effective phase stiffness" 
model <|35(1 • Precise description of the soliton structure is 
rather complicated and beyond the scope of this paper. 
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FIG. 9: Pinning of a Josephson vortex by a single pancake- 
stack row. 



III. PINNING OF JOSEPHSON VORTEX BY 
PANCAKE VORTICES 



with A' x (y) = dA x (y) jdy. For the maximum force we 
obtain 

F °^ 
xmax 47r 2 7 3 s 2 ln(A/ S /3 / ) 

with Cf — ma,Xy [A' x (y)] and (3f ~ 1. Numerical cal- 
culation gives Cf s» 1.4 and the maximum is located at 
Df~ 0.527s. The critical current which detaches JV from 
the row of pancake stacks with the period do is given by 

CVA 2 
7sa ln(A/s/?/J 

where jj = c$o/(87r 2 A 2 s) is the Josephson current. This 
expression is valid as long as the lattice period ao is much 
larger than the Josephson length 7s. Otherwise inter- 
action with several pancake rows has to be considered, 
which will be done in the next section. 



B. Pinning by dilute pancake lattice (A < ao < 7s) 



Dynamic properties of JVs can be probed either by 
applying transport current along the c direction or by 
studying ac susceptibility for magnetic field, polarized 
along the layers. Although there are numerous experi- 
mental indications that pancake vortices strongly impede 
motion of JVs^22i2&2i no quantitative study (theoretical 
or experimental) has been done yet. In this section we 
consider a pinning force, which is necessary to apply to 
the JV to detach it from the pancake vortex crystal. We 
consider the simplest case of a dilute pancake lattice, 
a > A, which allows us to neglect the influence of pan- 
cakes on the JV core, and small concentration of JVs, so 
that we can neglect influence of the JV lattice on the pan- 
cake crystal (e.g., formation of phase-separated states). 



A. Pinning by a single pancake-stack row (ao > 7s) 

We consider first the simplest case of an isolated 
pancake-stack row crossing JV (see Fig. |§J) and estimate 
the force necessary to detach JV from this stack, assum- 
ing that its position is fixed. The consideration is based 
on the crossing energy of JV and isolated pancake stack, 
calculated in Sec. Ill Bl Let us calculated first the force 
necessary to separate JV from an isolated pancake stack. 
Using Eq. (|17fl we obtain for the force acting on JV from 
the pancake stack located at distance y from the center 
of JV core 



When several pancake-stack rows fit inside the JV core 
(but still ao > A) interaction energy of JV per unit length 
with the pancake lattice, Eji(y), can be calculated as a 
sum of crossing energies (|17|) . 



Eji(y) 



— y^E x (y~ nb ) 



(67) 



47r 2 7 2 sa ln(7s/A) f 



y - mb c 

7s 



where bo = ^/3ao/2 is the distance between the PV rows. 
This expression has a logarithmic accuracy, i.e., we ne- 
glected a weak y-dependence under the logarithm and 
replaced the dimensionless function under the logarithm 
by its typical value. Using Fourier transform of E x (y), 
Exip) = J dyexp(-ipy)E x (y), we also represent Eji(y) 
as 



Eji(y) =n v ^ E x 



2ns 



exp % 



2irm 



-y 



In the case bo <C 7s we can keep only m = 0,±1 terms 
in this sum, 

Eji{y) w Ejt(0) + 2n v E x (^j cos (^^j 

Using interpolation formula (|21|l for A x (y), we obtain 
the asymptotics for bo < 7s 



F* = - Ty EAy) 



47r 2 7 3 s 2 ln(A/wi(y)) 



A' 



(65) 



Eji(y) - Eji(0) 



-0.35 



7111(73/ A) V bo 



(2<K V \ 

x cos — — exp 
V 60 J 



£ (68) 



-5.82 



7s 
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the relation between the JV-pancake interaction and 
disorder-induced pancake pinning. Moving JVs either 
can drag the pancake lattice or they can move through 
the static pancake lattice. Slow dragging of pancake 
stacks by JVs at small c axis fields have been observed 
experimentally. 17 As the JV-pancake interaction force de- 
cays exponentially at high c axis fields in the case A < 7s, 
one can expect that JVs will always move through the 
static pancake lattice at sufficiently high fields. 



FIG. 10: The field dependence of the dimensionless pinning 
force of JV by dilute pancake vortex crystal 



To compute the field dependence of the critical current, 
we represent the force acting on JV from the pancake 
lattice in the scaling form, 

F Jpi(y) ~ — f ( — , — 



(47r) 2 (7s) 4 ln(A/s) \7s'7s 



with 



F(y~b ) = J--^^A x (y-j6o). 
b dy ^ 

The critical pinning current is given by 

v^A 2 - 
Jjpl=3J 2( 1S )nn(X/sf c{bo) 

with 



(69) 



A. Dragging pancake lattice by Josephson vortices 

When moving JVs drag the pancake lattice, one can 
obtain simple universal formulas for the JV viscosity co- 
efficient and JV flux-flow resistivity. The effective vis- 
cosity coefficient per single JV is connected by a simple 
relation with the viscosity coefficient of pancake stack per 
unit length, t] v , 

B x r/j V = B z rj v . 

Therefore, the JV flux-flow resistivity, p^(B) = 
$ B x /(c?r)jv), is given by 



Pff 



'B) 



c 2 ri v B z 



(71) 



As a consequence, we also obtain a simple relation be- 
tween pj^(B) and p a ^(B z ) 

p} f (B) = P f f (B z )§. 



T c (b Q ) = max[J r (y, b )) 

y 

Numerically calculated dependence of T c versus reduced 
field (7s) 2 £/$ = \/3/26g. Maximum F cm3X » 1.15 is 
achieved at B « 0.26$ /(7s) 2 (bo = 2.1js). Therefore 
the maximum pinning current can be estimated as 



Jjp r 



A 2 



1 J J 



(7s) 2 ln(A/s) 



(70) 



For typical parameters of BSCCO this current is only 5- 
10 times smaller than the maximum Josephson current, 
i.e., it is actually rather large. 

An exponential decay of the pinning energy l|Q8[) and 
force holds until the row separation bo reaches A. At 
larger fields one have to take into account shrinking of 
the vortex core. Because the number of pancake rows in 
the core is almost constant at high field, the exponential 
decay will saturate at a finite value ~ exp(— Cjs/X). 



IV. PANCAKE VORTICES AND VISCOSITY OF 
JOSEPHSON VORTEX 

A sufficiently large c axis transport current will drive 
JVs. Two dynamic regimes are possible, depending on 



B. Josephson vortex moving through pancake 
lattice 



Consider slow JV motion through the static pancake 
lattice. JV motion along the y direction with velocity V y 
induces traveling pancake displacement field u n (y — V y i). 
For slow motion, u n (y) is just the static displacement 
field around JV. Contribution to the energy dissipation 
caused by these displacements is given by 



W « r] p n v 



Vp 



•E 



J dr (V y u n y 



where r\ v is the pancake viscosity coefficient. Therefore 
the JV viscosity per unit length is given by 



VJV = r\ v n 



V ^2 J dr(\7 y u n f 



(72) 



Using relation between the vortex phase and displace- 



ment field V 



y'vvn 



2nn v u n we obtain an estimate 



n v y j j dy{VyU n ) 



2.7X 3 y/E; 

(is) 3 



C 



3/2 



hi(r- cut /r u ,) / 
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with C = (cos(0„ + i — </>„)) < 1, and derive an approxi- 
mate formula for the JV viscosity coefficient 



Vjv ~Vv~ 



2.7A 3 

a (7s) 3 



C 



3/2 



In- 



(73) 



with r\ v = r/p/s is the viscosity of pancake stack per unit 
length. Because the JV viscosity rjjy is proportional to 
the pancake-stack viscosity r\ v , there is a relation between 
the flux- flow resistivity of JVs {p c jf) and the flux- flow 
resistivity of pancake vortices (pff) 



p)j_ _ pff a 0.37( 7 s) 3 ( Hr cut /r v 
B x ~ B z s A 3 I 



C 



3/2 



(74) 



If we use the Bardeen-Stephen formula for the in-plane 
flux-flow resistivity, pjj « PabB z / H C 2, and an estimate 
for the c axis flux-flow resistivity at B z = 0, Pf°f « 

(I67 3 s 2 B x /$o)p a b^ we can also obtain relation between 
p c ff (B z ) and pf f 



Ph^0.15pf f ^ 



C 



3/2 



(75) 



From this estimate we can see that at a ~ A the flux-flow 
resistivity for JVs slowly moving through the pancake 
lattice is about four order of magnitude (factor (£/A) 2 ) 
smaller than the flux-flow resistivity of free JVs. We see 
that even though the critical force becomes exponentially 



small at high fields, pancakes still very strongly hinder 
mobility of JVs. 
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APPENDIX A: MAGNETIC INTERACTION 
BETWEEN PANCAKE ROWS 

In this appendix we derive several useful representa- 
tion for the magnetic interaction between pancake rows 
in different layers. The interaction energy between two 
pancakes in Fourier and real-space is given by Eq. J7J). 
We will use also this interaction in the mixed represen- 
tation 



kl 



X 2 k 2 



The interaction energy between the pancake rows per 
unit length is given by 



UMr(x,y,n) = - y^U M (x - ma,y,n) 
a * — ' 

m 

a* J J 2ir a \ a 



From this equation we obtain the following integral representation for Umt(x, n ) 



UMr(x,y,n) 



seo 



S n - 777 exp I - 



2A 



s\n\ 



In 



2ttx 

1 — 2 cos exp 

a 



exp 



Air\y\ 



ex p i^fx + iJ\- 2 + (3si) usn - J \~ 2 + \/T+tf\y 



l + (¥) u a Wl 



This representation can be used to derive large-y asymptotics of Umt{x, y, n), 



U Mr (x,y,n) — 
a 



2irs 



exp % 



s\n\ 



S n exp — 2 cos exp 



2ttx 



2n\y\ 



A- 2 + 



1/4 
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aX\J A 2 + {2ir/a) 2 \y\ 3> 1. Therefore the interaction between the pancake rows decays exponentially at y, sn> A, a 
Interaction of pancake row with stack of rows is given by 



U Ms (x,y) = 2^U Mr (x,y,n) 

71=1 

ttJ 



■In 



a 



2irx ( 2th/ \ ( 4th/ \ 

1 — 2 cost ) cxp + cxp 

a \ a J \ a J 



+ 51 = 

a i V /(a/A) 2 + (2 7 r0 2 

In particular, the potential created by pancakes belonging to the same row stack ("cage potential") is given by 



t^cagc = U MS (X,0) 



Air J 



2 t 00 



V a 



2nl 



^(a/Xf + (2^) 2 



APPENDIX B: LOCAL CONTRIBUTION TO 
JOSEPHSON ENERGY DUE TO MISMATCH OF 
PANCAKE ROWS IN NEIGHBORING LAYERS 



We measure all distances in units of the lattice constant 
a. The x-averaged phase mismatch is given by 



We consider two pancake rows in neighboring layers 
with period a <C 7s shifted at distance v with respect 
to each other in the direction of row (x axis). In zero 
order with respect to the Josephson coupling these rows 
produce the phase mismatch ip v (x, y) between the layers, 



x— m + v/2 



ip v (x,y) = 'y~* j (arctan 

™ v y 



arctan 



x— m — v/2 

y 



<Pv(y) = ^sgn(y). 



The total phase approaches y> v {y) at y > l/2n. 

Separating ip v (x,y) from the total phase, we can rep- 
resent the Josephson energy as 



Tj = Eja 2 j d 2 r [1 - cos (ip v (x,y) + ip(x, y))] 

= E. r a 2 J d 2 r [1 - cos {Tp v {y) + <p(x, y))] + L x £ Josc (v, <p) 

where (f(x,y) is the smooth external phase and the local Josephson energy £j osc (v,ip) per unit length is defined as 

Eja 2 



■'Jose 



(v,<p) 



Eja f 

^ — / dxdy [cos (<p + <p v (x,y))- cos (<p + tp v (y))] 



(Bl) 



= 2Eja cos {if) J 



where 



J(v) 



1/2 poc 

dx dy [cos {ip v {y)) - cos {ip v {x,y))}. 
■1/2 Jo 



(B2) 



In £j sc(v, <f) we can neglect weak coordinate dependence of the external phase and replace it by a constant ip. At 
\v\ < 1/2 the ground state for fixed v corresponds to ip = 0, while at 1/2 < \v\ < 1 the ground state corresponds to 
ip = it. £j osc {v, ip) has a symmetry property £j osc {l — v, ip) = —£ Josc (v, if). The integral over y in J{v) is converges at 
y < l/2n. This allows us to consider a single row separately from other rows and neglect the coordinate dependence 
of the "external phase" ip. 
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Using the complex variable z = x + iy one can derive a useful expression for (p v (x, y): 



(p v (z) = Im I In JJ 

\ m 



z — m — v/2 
z — m + v/2 

Going back to the (x,y) representation, we obtain 



Im In 



sin (n(z — v/2)) 
sin (tt(z + v/2)) 



and 



. . tan (n(x + v/2)) tan (n(x — v/2)) 

tp v (x, y) — arctan ; arctan ; 

tanh Try tanh 7ry 

cosh 2ny cos ttv — cos 2?tx 



cos ip v (x,y) 



-^/(cosh 27ry — cos 27nc cos 7ri;) 2 — (sin27rxsin7n;) s 



Integral l|B2|) can now be represented as 



J(v) = 



1/2 



1/2 



dx I dy 



cos (nv) 



cosh 2ny cos ttv — cos 2itx 



y / (cosh27ry — cos 2irx cos 7rv) 2 — (sin27rasin7ri;) 2 



We obtain an approximate analytical result at small v, v <C 1. Simple expansion with respect to u produces loga- 
rithmically diverging integral. To handle this problem we introduce the intermediate scale yo, v <C yo <C 1, and split 
integral J into contribution coming from y > yo (jj>) and y < yo (J7<). At y > yo we use small- v expansion and 
obtain 



J> = M 

(ttvY 



dy 



1/2 



dx 



sinh 27T?/ 



(cosh 2iry — cos 2nx) 



2 

„2 



exp(— 27T?/) 



dy- 

yo sinh 2?ry 
1 



7TU 

w In . 

4 47ryo 

In region y < yo we can expand all trigonometric functions and obtain the integral 

4(x 2 + y 2 ) - v 2 



J< 



dy 



v 
~2 



,1/2 




/ dx 


1 - 


Jo 






dy / 


dx 


Jo 





a/(4 (a; 2 + y 2 ) + v 2 ) 2 - (4xv) 2 
1 



x 2 + y 2 - 1 



yj{x 2 + y 2 - l) 2 + Ay 



r,2 



I 

with y — 2y/v and x — 2x/v. Because we only interested We also calculated function J(v) numerically for the 
in the main logarithmic term, we can extend integration whole range < v < 1/2. The result is shown in Fig. 1111 
over x up to oo. The obtained integral can be evaluated and is described by approximate interpolation formula 



J<t * V fin £ + 1.58 
4 V v 

Adding J7> and J7<, we obtain 



J{v) 



1 — COS TTV 



47T 



(B5) 



J{v) 



-v In 



0.9f 



(B3) 



x In 



1 



1 — COS TTV 



0.379 costtw + 0.076 cos^ ttv 



and 



Ejosc(v, tp) Pa -Eja cos (p) v 2 In - - 0.95 (B4) 



This result was used in numerical calculations of the JV 
core structure. 
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FIG. 11: Dimensionless function J(v) which determines the local Josephson energy. Dashed line is small-u asymptotics 1B3II . 



APPENDIX C: CONTRIBUTION TO JV 
ENERGY COMING FROM LOCAL JOSEPHSON 
TERM AT 7 > A/s 



energy is given by 



In the limit of very weak coupling the correction to 
reduced JV energy l|57ll coming from the local Josephson 



6£ 



nV3 



JV 



VI + B Z /B X n i 



^ f° j~/ x2, 0-39 , , 

> / dy(v n+ i~v n ) In, r cos(0„ +x - 0„) 

tr'^J-oo \v n+1 -v n \ 



(CI) 



with y = y/Aj = yl + B z /B\y/js. We will focus only on the regime B z ^> £?a, where this correction can be 
noticeable. In this regime reduced row displacements are connected with phase gradient by relation 



Vn{y) 



27a7 v ^" 



and the correction reduces to 



5S.JV ~ — Bz [ BX ,2 y2 [ d V (Vj/^n+i - V a 0„) In 
8im v Hs) Z J K ' 



Vy^n+l — WyC, 



COS(0„+1 - 4> n ). 



8nn v (7s) 
Using numerical estimates 

00 roc . 2 

/ dy cos {<f) n+1 - 4> n ) \ Vy4> n+1 - V y (f} n ) ~ -2.4 

n=-oo"'- 00 

°° />oc 

^ dy COS (0 n + l - ) (Vy0 n+ 1 - Vy 



In 



1 



n— — 00 



4.7 



obtained with the JV phase 4> n (y) (fHE|) . we obtain 



/BT 0.38A 2 

OCJV ~ — 7T4 / — 5 m 

V Bz (7s) ln(a/r«,) 



0.1 7s / a 



(C2) 



(C3) 



As we can see, the correction is smaller than the reduced JV energy at B z ^> B\, Ejy w n-J 1 B\/B z ln(a/s), by the 
factor — A 2 /(7s) 2 . 
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